RU-98-28 



A practical implementation of the Overlap-Dirac operator. 



Herbert Neuberger 
neubergOphysics . rutgers . edu 



Department of Physics and Astronomy 
Rutgers University 
Piscataway, NJ 08855-0849 



Abstract 

A practical implementation of the Overlap-Dirac operator jg presented. The 

implementation exploits the sparseness of H and does not require full storage. A simple 
application to parity invariant three dimensional SU{2) gauge theory is carried out to 
establish that zero modes related to topology are exactly reproduced on the lattice. 
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The preservation of global chiral symmetries on the lattice has been a long standing 
problem. During the last five years a solution has been found but in the initial stages of 
the finding practical applicability seemed limited to two dimensions (2d). The overlap was 
developed starting from [1] , which in turn was motivated by [2] and [3] . The formulation 
of vector-like theories which is the main focus of this letter can be found in section 9 of 
[4] . Work on 2d examples was successful [5] and added confidence in the basic idea. This 
induced a subsequent search for simplifications. One year ago a compact action was found 
representing exactly massless fermions on the lattice [6]. This compact action is induced 
by integrating out all but one of the infinite number of fermions that are equivalent to 
the overlap. The fermions that are integrated out are heavy and the "miracle" is the 
relative simplicity of the action for the remaining massless fermion. That action confirms 
to a criterion introduced many years ago by Ginsparg and Wilson (GW) [7] to assure 
continuum-like Ward identities on the lattice. The compact action produced by the overlap 
is the first known acceptable solution to the GW requirement. To date, there are no other 
explicit solutions. Related work can be found in [8,9]. 

In spite of its compact form, the overlap-Dirac operator isn't easy to work with nu- 
merically and its practical usefulness in dimensions higher than two has been unclear until 
now. One purpose of this letter is to present a new procedure to simulate the overlap-Dirac 
operator on the computer. This procedure holds the promise to be practicable in 3 and 4 
dimensional gauge theories. 

One of the central issues in gauge theories is the spontaneous breakdown of global 
chiral symmetries. This issue could not be addressed in a clean way on the lattice until 
now, but the new development makes this possible. 

In three dimensions there are interesting models that are parity invariant and have 
analogues of four dimensional global chiral symmetries [10]. Intensive investigations have 
indicated that spontaneous breakdown of these symmetries takes place if the (even) number 
of fiavors is small enough. 

In the continuum spontaneous chiral symmetry breakdown is correlated with an en- 
hancement in the spectral density of the Dirac operator around zero. Clean checks of this 
were impossible on the lattice before the new actions. 

In four dimensions it has been suspected that configurations consisting of superposed 
instantons/anti-instantons are responsible for this accumulation of almost zero eigenvalues. 
In two dimensions it is quite likely that this is actually true [11]. In three dimensions there 
are no instantons, but there is another source of low lying levels. Level crossings at zero 
must occur in backgrounds connecting topologically distinct pure gauge configurations. 
These crossings also happen in four dimensions. Until now, the level crossings could not 
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be realized on the lattice since the needed symmetries were broken by the regularization. 

Using the new action and the new procedure I shall exhibit below a level crossing in 
three dimensional SU{2) gauge theory with two flavors. At the moment I see no major 
difficulties left to overcome on the way to a full dynamical simulation of this exactly 
"chirally symmetric" model. Generalizations to four dimensions also no longer appear 
prohibitive. It is possible that this new procedure will revolutionize the way fermions are 
treated on the lattice. 

The compact overlap-Dirac operator introduced in [6] is 

^^i±^, (1) 

where H — ^^Dy^. There is some freedom in choosing Dy^. The simplest choice is to take 
Dw as the Wilson-Dirac operator with hopping parameter k set at re = 2d-2 • This is the 
choice adopted henceforth. 

The first attempt to use D directly, rather than the original overlap formula [4], was 
made in two dimensions for a U{1) gauge theory by T-W Chiu [12]. He used a Newton 
iteration to find \fH^. His method required storage of the full matrix H in memory. It is 
easy to bypass the computation of V H'^ and establish an iteration for e{H) directly: 

^ Xk + ^]- k ^0,1,2,...; Xo = H. (2) 



Observe now that the iteration can be "solved" by the replacement Wk = i^xl • ^ assume 
that Xq has no eigenvalue equal to and if this is true it will be also true of all Xk- 
Equation (2) immediately leads to Wk = (Wq)^ . Setting n = 2^~^ ^ we write the solution 
to (2) as Xk = fn{H). The function fn{z) is given by 

f I.) - + (^^ 

•^"^ ^~ (l + ^)2n+(l-^)2n- 

fn{H) is a rational approximant to e{H) of the Pade type. The following identity is easily 
derived by calculating all the poles of /n(-z) and their residues: 



z 1 
n^2;2cos2^(s- |) + sm ^{s - ^ 



fn{H) is relatively easy to calculate using, for example, a CG algorithm for the inversions. 
Equation (4) has been previously derived by applied mathematicians [13] in a different 
context. 
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The computational cost is not much higher than a single inversion since the inversions 
for all s are related by shifts [14]. The convergence of the CG iteration is controlled by 
the s = 1 term. There is no reason to require n to be a power of 2 any more, and fn{H) 
is a truncation of e{H) for any n > 1. In practice, below, I used n = 48 to get ten digits 
precision when acting on random vectors. Essentially all the precision required of the CG 
was maintained. Since is odd all that one needs to check is how close ^^'^"[by^^^^ is to 
unity for randomly chosen vectors h. Since e{\H) = e{H) for positive A we also have scaled 
versions of /„, with slightly different convergence properties. Since physics is concentrated 
in the spectrum of H close to zero it is not recommended to choose a A to ensure optimal 
overall convergence. For simplicity, I have kept A = 1 in what follows. 

The formula for H appropriate for three dimensions was obtained by dimensional 
reduction along direction 4 from four dimensions [15]. The boundary conditions obeyed 
by the fermions are picked periodic in direction 4 and anti-periodic in directions 1,2,3. I 
used the real form of H introduced in [16]. Starting from a single Dirac fermion in four 
dimensions we end up with two Dirac fermions in three dimensions. Let us denote ^^e{H) 
by V. V is orthogonal, and obeys 75F75 = 74F74 = V'^ and consequently [7475, V] = 0. 
The first two identities are of GW type. The three properties together constitute the 
lattice realization of the global "chiral" SU{2) symmetry known in the continuum. The 
reality of V implies that parity is conserved at the action level. In the continuum, it is 
believed that parity does not break down spontaneously, but the global SU{2) symmetry 
does. The breaking is accompanied by two massless "pions" . 

The crossing we wish to see occurs when we interpolate smoothly unit link variables to 
a lattice pure gauge configuration of special type. The gauge transformation defining the 
latter is a discretized form of a known nontrivial continuum map from the three torus in 
the continuum to the SU{2) group manifold [17]. The interpolation parameter is denoted 
by t and goes from to 1. There is a built in symmetry in the configuration under t — > 1 — t. 
Thus, the crossing is expected exactly at t= .5. 

To see the crossing I compute the lowest few eigenvalues of 1 + ^ (F + V'^) by a CG 
based variational method [18]. The lowest state is found to be doubly degenerate. It 
corresponds to two conjugate eigenvalues of V, 6=*=*^. For t close to .5, 9 is close to tt, and 
the crossing can be easily seen from a plot of the square root of the eigenvalue 1 + cos d{t) 
as a function of t. The flow is shown in the Figure 1. 

The crossing takes place on lattices as small as 6^ and has been confirmed by the 
indirect method of [17]. That method is based on the Herzberg - Longuet-Higgins effect 
[19]. Without any effort I increased the volume to 16^. At this point His a 32, 768 x 32, 768 
real matrix. In QCD it has been claimed by Lepage and co-workers that physical results 
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can be obtained on lattices as small as 6^ [20]. Such a lattice has a complex fermionic 
matrix of size 15, 552 x 15, 552. This is comparable to our example. 

Many refinements of the procedure outlined here are possible. Lurking behind the 
scenes are Chebyshev polynomials and they have appeared in this context before, in the 
analysis of another truncation of the overlap, studied in [21,22]. The light fermions in the 
latter truncation are usually referred to as domain wall fermions. There are many reasons 
why I suspect that the method presented here is superior to domain wall fermions. Space 
does no permit a more thorough discussion at this time. 

In this letter a new and direct approach to treat the overlap action on 
potentially realistic lattices in three and four dimensions has been defined. The feasibility 
of a serious numerical study in d = 3 of a certain class of interesting models has been 
established. These results open up many directions for future research. 

Flow in smooth gauge field 
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Figure 1. The level crossing. 
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